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Abstract 

We have performed a deep optical imaging survey using a narrowband filter 
{NB921) centered at A = 9196 A together with i' and z' broadband filters cov- 
ering an 814 arcmin^ area of the Subaru Deep Field. We obtained a sample of 73 
strong A^i?921-excess objects based on the following two color criteria; z' — NB921 > 1 
and i' — z' > 1.3. We then obtained optical spectroscopy of nine objects in our 
iVi?921-excess sample, and identified at least two Lya emitters a.tz = 6.541 ±0.002 
and z = 6.578 ± 0.002, each of which shows the characteristic sharp cutoff together 
with the continuum depression at wavelengths shortward of the line peak. The latter 
object is more distant than HCM-6A a.t z = 6.56 and thus this is the most distant 
known object found so far. These new data allow us to estimate the first meaningful 
lower limit of the star formation rate density beyond redshift 6; Psfr ~ 5.2 x 1O~^M0 
yr~^ Mpc~'^. Since it is expected that the actual density is higher by a factor of 
several than this value, our new observation reveals that a moderately high level of 
star formation activity already occurred at z ^ 6.6. 

Key words: cosmology: observations — early universe — galaxies: starburst — 
galaxies: formation — galaxies: evolution 

1. INTRODUCTION 

Probing the star formation activity in galactic or subgalactic systems in early universe is 
important for understanding both the history of galaxies and the origin of cosmic reionization 
(e.g., Loeb, Barkana 2001). Recent advance in deep optical imaging capability with 8-10 m 
class telescopes enabled new searches for star- forming galaxies beyond redshift 5. In particular, 
imaging surveys using narrow-passband filters have proved to be an efficient way to find such 
galaxies (Ajiki et al. 2002; Cowie, Hu 1998, Hu et al. 1996, 1999, 2002; Kudritzki et al. 2000; 
Rhoads et al. 2001; Steidel et al. 2000; Taniguchi et al. 2003). Indeed, the most distant Lya 
emitter (LAE) known to date, HCM 6A at z = 6.56 was discovered by using this technique 
(Hu et al. 2002). However, surveys for emission-line galaxies with narrow-band filters have 
an intrinsic limitation in redshift coverage and hence the survey volumes are not often large 
enough to ensure robust for success. In order to increase survey volumes and to reach faint 
limiting magnitude, we need wide-field CCD cameras on 8-10 m class telescopes. Suprime-Cam 

* Based on data collected at the Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 

This work has been done with a collaboration between the Subaru Deep Field Project led by the association 
of builders of the Subaru telescope and the common-use, Intensive Program (S02A-IP-2) led by Y. Taniguchi. 
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(Miyazaki et al. 2002) mounted at the prime focus of the 8.2 m Subaru Telescope (Kaifu 1998) 
on Mauna Kea, Hawaii, provides a unique opportunity for wide- field (a 34' x 27' field of view), 
narrowband imaging surveys for emission-line objects at high redshift. 

The most distant LAE, HCM 6A aX z = 6.56 (Hu et al. 2002), is gravitationally amplified 
by a factor of 4.5 by the foreground cluster of galaxies Abell 370 at z = 0.37. Although the 
help of any gravitational lensing is highly useful in investigating faint high-z objects (Ellis et 
al. 2001; Hu et al. 2002), it is also important to search for high- 2; LAEs in a so-called blank 
field for an unbiased study. In an attempt to find star-forming objects at 2; ~ 6.6 in such a 
blank field, we have carried out a very deep optical imaging survey in the Subaru Deep Field 
(SDF) centered at a(J2000) = 13^^ 24°^ 21.4^ and 5(J2000) = +27° 29' 23" (e.g., Maihara et al. 
2001; Ouchi et al. 2003; Kashikawa et al. 2003). In this Letter, we report on our discovery of 
two Lya emitters at 2 ~ 6.5 - 6.6. 

2. OBSERVATIONS 

2.1. Optical Imaging 

In this survey, we used the narrow-passband filter, NB921, centered on Ac = 9196 A with 
a passband of AA(FWHM) = 132 A, corresponding to a redshift range between 6.508 and 6.617 
for the Lya emission. Optical imaging was made in the i\ z' , and A^B921 bands on a central 
34' X 27' area of the SDF with Suprime-Cam. Our direct imaging was obtained during several 
observing runs between 2001 April and 2002 May. The total integrations time was 4.7 hr, 
5.8 hr, and 5.0 hr for i', 2;', and A^i?921, respectively. The data reduction procedures are the 
same as those given in Yagi et al. (2002). The PSF FWHM of the final images is 0."90. Source 
detection and photometry are performed using SExtractor version 2.1.6 (Bertin, Arnouts 1996). 
The limiting magnitude (AB) for a ba detection on a l."8 diameter aperture is 26.9, 26.1, and 
25.7 for z', z', and A^i?921, respectively. 

For each object detected in the A^i?921 image, i', z' , and A^i?921 magnitudes are mea- 
sured on a common aperture of l."8 diameter. In total, 50,449 objects are detected down to 
A^i?921 = 25.7 (the 5o" limiting magnitude). The effective area used to search for A^i?921-excess 
objects is 814.3 arcmin^. The FWHM half-power points of the filter correspond to a co-moving 
depth along the line of sight of 40.9 h^^j Mpc. Thus a total volume of 202,000 h^^ Mpc^ is 
probed in our NB921 image. 

We select candidates of 2 ~ 6.6 LAEs imposing two criteria, z' — NB921 > 1 and i' — z' > 
1.3, on the above objects. There are 404 objects satisfying the first criterion alone. The latter 
criterion is used to reduce contamination from foreground objects. Objects at z ^ 6.6 have 
sharp breaks because of the strong Lya absorption at this redshift (Songaila, Cowie 2002) and 
are expected to exceed the adopted i-z criterion while low redshift galaxies invariably do not. 
This latter criterion is applied only to objects brighter than i' = 28.0 (^2(T limiting magnitude), 
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and we retain all the objects (satisfying z' — NB921 > 1) with i' > 28.0 in order not to miss 
possible faint LAEs. These selection procedures yield eventually a photometric sample of 73 
LAE candidates. Note that among the objects with i' > 28, those brighter than z' = 26.7 
(3(7 limiting magnitude) satisfy automatically i' — z' > 1.3, but we cannot obtain a meaningful 
constraint on i' — z' for those with z' > 26.7. 

2.2. Optical Spectroscopy 

In order to investigate the nature of LAE candidates found in our optical imaging survey, 
we obtained optical spectroscopy of nine objects in our LAE candidate sample using the Subaru 
Faint Object Camera And Spectrograph (FOCAS: Kashikawa et al. 2002) on 2002 June 7-9 
under a 0."55 - l."10 seeing condition. This spectroscopic sample contains the brightest three 
objects with A^i?921 < 24. However, the remaining six objects are randomly selected from the 
photometric sample and they have NB921 ^ 25 - 25.5. 

Our optical spectroscopy was made with the 300 lines mm~^ grating and an 058 order 
cut filter. The wavelength coverage was between 6000 A to 10000 A with a pixel resolution of 
1.34 A. The use of an 0.8 arcsec-wide slit gave a spectroscopic resolution of 9.0 A at 9200 A (or 
R ~ 1020). The spatial resolution was 0."3 pixel^^ by 3-pixel, on-chip binning. Spectroscopy 
of the brightest object in our LAE sample {NB921 = 23.01) was obtained in the long-slit 
spectroscopy mode. The exposure time was 1800 seconds. This source was quickly identified 
as an [O lll]A5007 emitter a.t z ^ 0.84. Other eight objects were observed in the multi-object 
spectroscopy (MOS) mode. We chose two fields that contained as many LAE candidates as 
possible; hereafter field 1 and field 2. We succeeded in obtaining spectra of three targets in 
Field 1 and five in Field 2. The same grating, filter, and slit as those in the long-slit mode were 
used in this MOS mode. We obtained twelve and six 1800s exposures for Field 1 and Field 2, 
respectively. We also obtained the spectrum of a standard star Hz 44 for flux calibration. 

3. RESULTS 

Our optical spectroscopy (see Figure 1) indicates that at least two objects are well- 
defined LAEs between z = 6.5 and z = 6.6 because their emission-line shapes show the sharp 
cutoff on the UV side together with the continuum depression at wavelengths shortward of 
the line peak. Two more objects also show the sharp cutoff at wavelengths shortward of 
the line peak. However, their continuum magnitudes are so faint that we cannot see firm 
evidence for the continuum break. Although they are probable candidates of LAEs at z ^6.5 
- 6.6, we do not include them in the later discussion^. In the spectra of the well-defined two 

^ Spectroscopic properties of other five objects are as follows. Three objects show a single emission line around 
A9200 A in our optical spectra. Since they show an almost symmetrical emission-line profile, we cannot 
identify them as LAEs solely basing upon our optical spectroscopy; they may be either an [O ii]A3727 source 
at z w 1.46 or a LAE at z « 6.6 (e.g., Stern et al. 2000). It seems also worthwhile noting that a symmetric 
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LAEs, SDF J132415.7+273058 in Field 2 and SDF J132418.3+271455 in Field 1 [panels (a) 
and (b), respectively], continuum emission appears to be present at wavelengths longer than 
the Lya peak. Combining these two spectra and applying a 3-pixel smoothing, we obtain 
the average spectrum in panel (c). The average continuum flux density between 1200 A and 
1210 A is (1.7 ±4.1) X 10~^9 ergs s'^ cm'^ A'^ while that between 1219 A and 1226 A is 
(6.4±3.7) X 10"^^ ergs s~^ cm"^ A^^. This difference is regarded as evidence for the continuum 
break in the spectra of these two LAEs. 

The redshifts of the two LAEs are estimated from the peak of the Lya emission 
line and the results are given in Table 1 together with line widths. It is noted that SDF 
J132418.3+271455 {z = 6.578) is more distant than HCM 6A at z = 6.56 (Hu et al. 2002) 
because error of our redshift measurement is ± 0.002. Therefore, SDF J132418. 3+271455 is the 
most distant LAE known to date. Thumbnail images of the two objects are shown in Figure 2. 
The NB921 images reveal that only SDF J132415. 7+273058 is spatially extended; its angular 
diameter (FWHM) is estimated as 1.2 arcsec. Correcting for the seeing spread (0.9 arcsec), we 
obtain an angular diameter of 0.8 arcsec, corresponding to 4.4 Hq ]. kpc at z = 6.541. 

4. DISCUSSION 

The SDF is a so-called blank field and there is no apparent cluster of galaxies known to 
date at low and intermediate redshifts in our field. Since, further, the lensing effect is expected 
to be small in this field, our survey results allow us to perform a simple statistical analysis of 
star formation activity in the investigated volume even though our sample is not large. 

First, we estimate the star formation rate of the LAEs at 2; ~ 6.6 by using the relation 
SFR{Lya) = 9.1 x 10-^=^L(Lya) MqJt"^ (Kennicutt 1998; Brocklehurst 1971). The observed 
Lya flux, the Lya luminosity, and the star formation rate, SFR(Lya), of each LAE are summa- 
rized in Table 1 where a fiat universe with ^matter = 0.3, = 0.7, and h = 0.7 with h = Ho/ (100 
km s~^ Mpc~^) is adopted. The average SFR obtained for the two LAEs is 7.1 ± 2.0 h^^ Mq 
yr~^, being comparable to those of LAEs at 2; ~ 5.1 - 5.8 (e.g., Ajiki et al. 2002). It should 
be mentioned that the SFRs estimated above are lower limits because it is quite likely that 
a blue half or more of the Lya emission may be absorbed by H I gas and dust grains in the 
galaxy itself and by the intergalactic H I gas (Miralda-Escude 1998; Miralda-Escude, Rees 
1998; Cen, McDonald 2002). The SFR based on the Lya luminosity tends to be underesti- 
mated by a few times or more than that based on the UV luminosity (see also Hu et al. 2002). 
Indeed, using the average UV continuum flux density between 1219 A and 1226 A for the 
combined spectrum shown in panel (c) of Figure 1 together with the relation (Kennicutt 1998), 
SFR{\JN) = 1.4 X 10~^^L^ Mq yr~^, where L^, is in units of ergs s~^ Hz~^, we obtain an average 

emission- line profile does not necessarily rule out the case of Lya emission (e.g., ESO 350-IG038 in Kunth 
et al. 1998). The remaining two sources are confirmed to be [O iii]A5007 emitters at z w 0.84. 
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value of SFR{VV) ~ 22/i(^^ M© yr'^ for the two LAEs. As for SDF J132415. 7+273058, we 
obtained J-band imaging using the InfraRed Camera and Spectrograph (IRCS: Kobayashi et 
ah 2000) on the Subaru Telescope on 2002 July 15. The integration time was 6480 seconds (the 
detail will be given elsewhere). The J magnitude (AB) is estimated as ~ 24.9. This photometry 
allows us to estimate the star formation rate at Arest = 1650 A for SDF J132415. 7+273058; 
SFR{l]V) ~ 36hoj Mq Yr~\ which is even higher by a factor of four than SFR[Lja). This is 
suggestive of dust obscuration at the bluest wavelengths. 

Now we can estimate the total star formation rate of 73 LAEs in our photometric sample 
using the equivalent width of NB921 flux. We have identified one LAE (SDF J132415. 7+273058) 
in the three brightest candidates and the other (SDF J132418. 3+271455) in the faint sample. 
Because of small number statistics, it seems modest to assume that approximately 22% (=2/9) 
of 73 LAE candidates are real LAEs at 2; ~ 6.5 - 6.6 in any magnitude range; note that 95 % 
confidence level based on the random sampling hypothesis ranges 8 % - 49 %. If we assume 
that all the 73 LAE candidates are true LAEs at 2; pa 6.5 - 6.6, we obtain nominally a total 
star formation rate of SFR^^^l^^^ = ATbh^^ Mq yr~^. Provided that approximately 22% of 
73 LAE candidates are real LAEs at 2; ~ 6.5 - 6.6, we can estimate the total star formation 
rate, ^Fi?totai ^ 0.22 x ^Fi?"™*"^ - Mq jt-\ Given the survey volume, 202,000 hoj 

Mpc^, we thus obtain a star formation rate density of Psfr — 5.2 x 10~^/io.7 Mq yr~^ Mpc~^. It 
should be reminded here again that we apply neither any reddening correction nor integration 
by assuming a certain luminosity function for LAEs. Further, we note that there are two 
more probable LAE candidates in our spectroscopic sample. Therefore, our estimate should be 
regarded as a robust and first meaningful lower limit for the star formation rate density beyond 
z = 6. Let us compare this value with previous estimates in Figure 3; note that we convert 
all the previous estimates to those in the cosmology adopted in this Letter. In conclusion, the 
present study shows unambiguously that the moderate star formation activity already occurred 
in the early universe beyond z = 6. 

We would like to thank the Subaru Telescope staff for their invaluable assistance. We 
would also like to thank Prof. Hy Spinrad both for useful comments for the earlier version 
of this paper and for encouragement. We also thank the referee, Daniel Stern, and another 
anonymous referee for their useful comments and suggestions. 
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Table 1. Properties of the two LAEs 



No. 


Name^ 


Optical AB Magnitude 
i' z' NB921 


z 


FWRM*" 
A km/s 


A 


F(Lya) 
10^^^ crgs/s/cm^ 


L{Lya) 
10*^ crgs/s 


SFR{Lya) 

Me yr"^ 


1 

2 


SDF J132415. 7+273058 
SDF J132418. 3+271455 


27.48 25.82 
28.52 26.66 


23.99 
24.98 


6.541 
6.578 


10.9 357 
< 9.0 < 290 


160t?f 
330-200 


2.06 + 0.18 
1.13 + 0.05 


10.02 + 0.09 
5.58 + 0.25 


9.1 + 0.8 
5.1 + 0.2 



^ The sky position, q:(J2000) and (5{J2000), is given in the name. 



^ Full width at half maximum of the Lyct emission line. 

^ Equivalent width of the Lya emission line at observed frame. 
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Fig. 1. The rest-frame UV spectra of two LAE candidates between 1200 A and 1230 A. Ttie spectra 
of SDF J132415. 7+273058 and SDF J132418.3+271455 are shown in panel (a) and (b), respectively. The 
combined spectrum of these two LAEs is shown in panel (c). A trial of the profile fitting with a combination 
between emission and absorption is shown by red line. 
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Fig. 2. Thumbnail images of the two well-defined LAEs, (a) SDF J132415. 7-^273058 and (b) SDF 
J132418. 3+271455. The size of each image is 10" x 10", and north is up and east is left. 
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Fig. 3. The star formation rate density (psfr) as a function of redshift z. Our estimate at z w 6.6 (large 
filled diamond) is shown together with the results of previous Lya searches at z ~ 3 - 5 (CH98 = Cowie 
& Hu 1998, KOO = Kudritzki et al. 2000, F03 = Fujita et al. 2003, and O03 = Ouchi et al. 2003). 
The previous investigations are shown by filled triangle (Gallego et al. 1996), open triangle (Treyer et 
al. 1998), open circle (Tresse & Maddox 1998), stars (Lilly et al. 1996), open pentagons (Hammer et 
al. 1999), filled squares (Connolly et al. 1997), filled circles (Madau 1998), and open squares (Pettini et 
al. 1999). Other results for Lyman break galaxies between z = 3 - 4 are also shown by filled pentagons 
(Steidel et al. 1999). 
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